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Abstract. Since astronauts on future deep space missions will be exposed to dangerous radiations, there is a need to accurately 
model the transport of radiation through shielding materials and to estimate the received radiation dose. In response to this need a 
three dimensional deterministic code for space radiation transport is now under development. The new code, GRNTRN, is based on 
a Green's function solution of the Boltzmann transport equation that is constructed in the form of a Neumann series. Analytical 
approximations will be obtained for the first three terms of the Neumann series and the remainder will be estimated by a rton- 
perturbative technique. This work discusses progress made to date and exhibits some computations based on the first two 
Neumann series terms. 

1.0 INTRODUCTION 

A recent National Research Council Report 
on the management of space radiation risk 
[1] highlights the need for an accurate and 
efficient three dimensional radiation 
transport code to determine and verify 
shielding requirements. According to this 
report, predictions derived from radiation 
transport calculations need to be tested 
using a common code for laboratory and 
space measurements that have been 
validated with accelerator results. Studies 
by Wilson et al. [2,3] have identified Green’s 
function techniques as the likely means of 
generating efficient high charge and energy 
(HZE) shielding codes that are suitable for 
space engineering and are capable of being 
validated in laboratory experiments. In 
consequence, a laboratory code designed 
to simulate the transport of heavy ions 
through one or more layers of material was 
developed at NASA Langley Research 
Center [2, 4, 5, 6, 7], It was based on a 
Green’s function model as a perturbation 
series with non-perturbative corrections. 

This early code used a scale factor to 
equate range-energy relations of one 
material thickness into an equivalent 
amount of another material, and proceeded 
to perform the transport calculations in the 
new material [8], This method proved to be 
acceptable for use with low-resolution 
detectors [6,9], but is unsuited for high- = Fj(s. b ,£i,E), (1) 

resolution measurements. Range and 
energy straggling, multiple Coulomb 


scattering and energy downshift and 
dispersion associated with nuclear events 
were also lacking from the prior solutions. In 
recent publications [10,11], it has been 
shown how these effects can be 
incorporated into the multiple fragmentation 
perturbation series leading to the 
development of a new Green’s function 
code GRNTRN (a GReeN’s function code 
for ion beam TRaNsport), GRNTRN has 
accurately modeled the transport of ion 
beams through multilayer materials 
[9,12,13] and unlike earlier codes it does not 
make use of range scaling. It is, however, 
deficient in that no account is taken of the 
variation of particle flux with angle since the 
code is purely one dimensional. The present 
paper strives to remove this deficiency by 
using generalizations of previous work to 
develop a fully three dimensional GRNTRN 
code and reports on progress made towards 
that end. 

2.0 BODY 

Consideration is given to the transport of 
high charge and energy ions through a 
three-dimensional convex region V , that is 
bounded by a smooth surface dv , and is 
filled with a target material. It is assumed 
that dV is subject to a boundary condition 
of the form 
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where FAx b ,SL,E) is a prescribed function, where 


and x b is a point on the boundary. It is 
required that ft-n(jCj,)<0 , where n(.v 6 ) is 
the unit outward normal at x 6 e dV , and the 
index j takes on values from 1 to N where 
jV is the number of ions in the model. 
^(x,ft,£) is the flux of j-type ions with 

atomic mass A t , at x e F t moving in the 

direction ft with energy E in units of 
MeV/amu. 

2.1 Transport Theory 

According to Wilson [14], the flux is given by 
the transport integral equation 
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where /? = x ■ ft , x'-x~{p~p')i 2 is the 
point where the ray through x in the 
direction ft enters V , SAE) is the stopping 

power, P t {E) is the nuclear attenuation 
function and <r k (£l,£l\E,E') is the 

production cross section for j-type ions with 
energy E and direction ft by the collision 
of k-type ions with energy £' and direction 
ft'. In addition E and E” are defined by 

Ej = E, (p - p \ E) ^ R ; 1 [Rj (E) + p - p '| , and 
E"=E j (p-p"-E) where R t (E) is the usual 
range-energy relation. 

The production cross sections used in this 
paper are based on S. R. Blattnig’s model, 
which is fully described in reference [17j. 
They are given by the approximation 
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and a jk {E k ) is the total cross section. In 
these equations, m p is the proton rest 

mass, y L and fi L are parameters associated 
with the Lorentz transformation from the 
fragment reference frame to the lab frame, 

<9 = cos '(ft -ft*) is the lab frame scattering 

angle, E s is the lab frame energy shift, and 
a the corresponding fragment momentum 
width [15,16], 

By introducing the field vector 

<D(x.ft./i) = [^(x.ft,£)j, (6) 

the fragmentation operator E 

|E-q>|,(x,ft : £) = XL> A 'j ( / /n ’ 

k>j 

■cr ; ,(ft.ft' ! £\£V,(x.ft , ,£’), (7) 


and the linear transport operator L 


PAE")SAE”) 

=— 4±—!-fm,E''), 

Pj(E) Sj(E ) 

the transport integral equation can be 
expressed in the operator form 


( 8 ) 


d> = G " F + Q L E <D 


( 9 ) 
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where Q represents the integral with 
respect to p" . Since Eq (9) is a Volterra 



type integral equation, it admits the in [10] and [11], it can be shown that the 

Neumann series solution [10] zero order Green's function takes the form 


<U = £(Q-L-H)" -C"-F=£C"-F (10) 

!P=0 7J-0 

where F is the boundary flux vector and for 

n> 1 , 

G” = (Q L 3) G" 1 . (11) 
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In this solution, the term G°-F represents 
the primary flux vector and the term G" -F 
represents the flux of the n ih generation of 
secondary ions produced. 

When the boundary condition (1) takes the 
special form 

2tt 

■5{E-E a )S(x b -x a ) (12) 

where 5 is the ‘surface delta function’ on 
dV , the solution of Eq. (2) is called the 
Green's function and is denoted by the 
symbol G jm (\,x 0 ,Sl,Sl 0 ,E,E 0 ) . Once the 

Green's function is known, the solution for 
an arbitrary boundary condition (1) can be 
obtained from the formula 

^(x,ti.£) = [G-Fl,(x ; a£) 

<? A {x,x 0 ,n 1 n 0 ,£,^)F i (x 0 ,n 0 A> ■ O 3 ) 

The summation is taken over k > j (instead 
of k > j) to account for the primary ion 
spectrum. 

2.2 The Zero Order Green’s Function 

The zero order Green’s function is the first 
term in the Neumann series (10) with the 
unit boundary condition (12). On taking 
account of energy straggling, as described 


where E m = E m (p -p\ E) and, by definition, 
E^EJp-p^^R^lRJEJ-ip-p^ 
is the mean energy at depth (p -p') g/cm 2 
of an m-type ion that entered the transport 
material with energy E 0 MeV/amu, and 

s m (p-p\E 0 ) is the corresponding energy 
straggling width. 

When the boundary condition takes the 
more general form (1), the primary flux is 
given by 

^;(v,fi,A') = [G 0 -F],(x,fi./,’) 

•Gj (*, x (1 . £2, P- . E, E 0 )F, . E 0 ), ( 1 5) 

which, in general, needs to be evaluated 
numerically. In the accelerator beam model 
described below however, Eq, (15) can be 
approximated analytically and a closed form 
expression obtained for the zero order 
primary flux. The result obtained in this case 
is called the broad zero order Green’s 
function. 

Since ion beam experiments play an 
important role in analyzing the shielding 
requirements against dangerous space 
radiations, there is interest in modeling the 
propagation of linear accelerator beams 
through potential shielding materials. A 
simple model can be constructed by 
assuming that the accelerator beam 
consists of m-type ions with mean energy 
E 0 MeV/amu and mean direction £2 0 . It is 
further assumed that the beam has 
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Gaussian profiles in both angle and energy 
and that it enters the material at points that 
are distributed in a Gaussian manner about 
the mean point of entry x 0 . In order to 
accomplish this, it is assumed that the 
boundary dV is defined by the single- 
valued, continuously differentiable 
parametric equations 

ciV = {x : x = x(u, v), u s <ti< w /j v s < v < v,.} , in 

which case the element of surface area is 
given by dS = |t)„x xd^dudv and the 
surface delta function is given by 

S < (x - x 0 ) = |9„x 0 x9„x 0 |”' S{u - u 0 )S(v - v 0 ) . 

The boundary condition (1) may then be 
assumed to take the Gaussian form 
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where x h = x b (u b ,v b ) , s x ,s a and s /; are the 
spreads in space, angle and energy 
respectively, and K a , K x are normalization 
constants. It should be observed that in the 
limit as s x ,s a ,s E — > 0 , the boundary 
condition (16) reduces to the Green’s 
function boundary condition (12). 

Equation (16) may now be substituted into 
Eq. (15) and the resulting integrals 
approximated by the mean value theorem 
and saddle point techniques discussed in 
[10]. The primary flux, which in this case is 
called the broad zero order Green's function 
G* m (x,x 0; ft.ft 0: £\.E 0 ) > is then given by the 
expression 




f r. 
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2.3 The First Order Green's Function 

The first order Green's function is given by 
the second term of the Neumann series (10) 
with the unit boundary condition (12). Since 
boundary condition (12) is a special case of 
the boundary condition (16), only the latter 
will be discussed. It may be recalled that 
the first generation fragment flux is 
determined by the formula 
G'=(Q L -S) G°. Therefore, on replacing 
G" by the broad zero order Green's 
function G* and expanding the result, it is 
found that the broad first order Green’s 
function is given by the expression 


G ! m (x. x,-,, ft. ft,,. E, E 0 ) = 



PjWSjiE) 



| 5 r , x i va v x i'|VT/,!,(A P\ \e q ) 


343 



' 1 2.v 


|g,-E,(A"-/?, ',£,)] 

T ^ ^ t r vl 


2 ^(A-AW 


(l-n.-o,,) 2 


-cxp<- 


where /?," = x"-£J 1 , x.'-x "-(/?, "-p.’JQ, is 
the point where the ray through x" , in the 
direction 12, enters v, E m =EJp, n - p,’,E,), 

K = KiP\ "“A \E 0 ) is the mean energy at 
depth (p p ] ’) g/cm 2 of an m-type ion that 
entered the transport material with energy 
E 0 MeV/amu, and $*(/>, "-p ) \E 0 ) is the 
corresponding energy width. 

The expression in Eq. (19) can be evaluated 
by numerical quadrature, but this is 
computationally expensive and therefore it 
is desirable to construct an analytical 
approximation. This can be done by making 
use of Taylor's theorem, the mean value 
theorem, and saddle point techniques as 
described in [10], and yields the result 
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p* is the root of the equation g Jm (p”) = 0 , 
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3.0 DISCUSSION 

To illustrate some of the theory presented 
above, the case in which the boundary flux 
consists of the m th component of the 
Galactic Cosmic Ray (GCR) spectrum is 
discussed. Since the GCR is isotropic and 
spatially uniform, boundary condition (1) 
takes the special form 

,E) = S M F m (E), (25 

where the spectra FJE) are broad 
functions of energy. These have been 
modeled by Badhwar and O’Neil (1995) 
and made available in tabular form at a 
number of solar maxima and minima 
between the years 1958 and 1989. 

In these circumstances, Eq. (15) for the 
primary flux takes the form 
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and, with the help of Taylor's theorem and 
the mean value theorem, may be further 
approximated as 


H. E) ~ S im H \ -£2 n(x ') | 
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The first generation fragment flux also 
simplifies and is approximated by the 
expression 

^ , (x.a/t) = |G , F|,(x.a£) 

=// |-£2 n(x’)| 

\~G)Ap,p\E,E 0 )F m {E 0 )dE 0 (28) 


where G'f m (p,p',E,E 0 ) is the special case of 
G%(p,p\E,E 0 ) for which s E =0. 

3.1 GCR on a Half-Space 

In the first example to be considered, the 
target is an aluminum solid that occupies 
the half-space V = {(x,y,z):z>0} whose 
boundary dV is the xy- plane . The 
Cartesian components of the vector £2 are 
given by £2 = (sin^cosor.sin^sinor r cosy) and 
are therefore completely determined by the 
polar angles y and a . The measured GCR 

56 Fe flux associated with the 1977 solar 
minimum [18] provides the boundary 
condition (25). The primary flux (27) then 
takes the form 


f J (*,a,E)~S Jm H[cosy\ 

P m (E m )S m (E m ) F E 
PJE) SJE) * 


(29) 


where E m =E m (z/ cos y,E) and 0<y</r/2. 


Since the field is axisymmetric about any 
line parallel to the z-axis , the fluxes of 
interest are functions of z,y, and E only. 
Figures 1 and 2 show the variation of the 
56 Fe primary flux with y and E at the 
points (0,0,0) and (0,0,5) respectively. 

Figure 3 provides a similar illustration for the 
first generation ' 6 () fragments at the point 
(0,0,5). 



Fig 2. 



Fig 3. 
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3.2 GCR on a Circular Cylinder 

In the second example to be considered, 
the target is the solid aluminum cylinder 
x 2 +y 2 < 16- , 0 < 2 < 36 and again the 
boundary condition is provided by the 
measured GCR 56 Fe flux associated with 
the 1977 solar minimum. 

Figures 4 and 5 show the variation of the 
56 Fe primary flux with y and E at the 
points (0,0,0) and (0,0,18) respectively, 
where the field is axisymmetric. Similar 
results for the first generation U, G fragment 
flux are shown in Figs. 6 and 7. 

The remaining figures deal with the flux at 
the point (14,0,9) where the field is no 
longer axisymmetric. Figures 8 and 9 show 
how the 56 Fe primary flux varies with y and 
E when or = 0° and or = 90° respectively. 
Figures 10 and 11 provide a similar 
illustration of the ]6 () fragment flux. 
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FIG 6. 


FIG 7. 


FIG 8. 


FIG 9. 


FIG 4. 


FIG 5. 
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FIG 10. 



FIG 11. 

4.0 CONCLUSIONS) 

In this work, some progress toward the 
development of a fully three dimensional 
deterministic code for space radiation 
transport was discussed. Approximations 
were obtained for the first two terms of the 
Neumann series solution of the transport 
integral equation. The results were then 
illustrated by exhibiting the primary flux and 
the first generation 16 0 fragment flux in a 
half-space, and in a circular cylinder, due to 
the 5S Fe component of the GCR. 

In future work an approximation for the third 
Neumann series term will be obtained and 
the Neumann series remainder estimated 
by a non-perturbative technique. Predictions 
made by the code will then be compared 
with the results of laboratory beam 
experiments and measured space data. 
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I n t rod ut t ion / M otivation 

Radiation in Space 


3 Main Types of Radiation: 

® Galactic Cosmic Rays (GCR) 

3 energetic charged particles 
3 penetrating power 

o Solar Particle Events (SPE) 

3 energetic protons and alpha 
particles 

® not likely to fragment 

o Particles Trapped in Radiation 
Belts 



Figure: Earth's magnetosphere 
and its interaction with the sun. 
( NASA ) 
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I n t rod ut t ion / M otivation 

Damage from Radiation 



Figure: DMA Damage due to heavy ions.^ nasa ) 


Consequently, steps must be taken to ensure astronaut safety by providing 
adequate shielding; 
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Transpoft Theorv 


Future Space Programs 



Figure: Solar Particle Events. 

( NASA ) 


® Require accurate and efficient 
methods for radiation transport - 
9 Determine and verify shielding 
requirements 

a National Research Council Report - 
9 Predictions need to be validated 
9 Use a common code for lab and 
space measurements 
9 Capable of being validated with 
accelerator results 

a Green's function techniques - 

9 Likely means for space engineering 
and lab experiments 
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Transport Theorv 


Volterra Integral Equation 



where 


n = 

Direction of propagation 

~Sj(E) = 

Energy lost per unit path length per unit mass 

aj(E) = 

Macroscopic absorption cross section 

*(X, ff.E'H 

Flux of k-type ions 


Double d if fe re nt ia 1 prod u ct io n c ross sect ion 

Pj(E) - 

Nuclear survival probability 
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Transport Theorv 

Green’s Function 



Figure: Flux of particles through a material. 


[ The Solution: Green’s Function 

1 

E()\ 

(2) 
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Transport Theorv 


Neumann Series Solution 




Figure: Atmospheric air showery Pierre Auger Observatory Team ) 
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Zero Order Green's Function 

Primary Flux 


[ The primary flux can 

be obtained by 


<p](x.n,E) = 

[G°F ]j(x,Q,E) 


= 

t dx o / dQ o / dE 0 
JdV J 4tt Je 



■Gj(x, xo, fi. O 0 , E, E 0 )Fj(x 0 . Qq. Eo), 

(4) 

- 


which needs to be evaluated numerically for some situations. 
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First Order Green's Function 

First Order Green’s Function 


= [G° + G 1 + G 2 + G 3 + , ] F 



Figure: Unit Boundary Condition. 
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Results 


Galactic Cosmic Rays (GCR) 


Figure: Galactic Cosmic Ray Distribution. (^5.^) 



^/(X|>i F) — 

where the spectra F m (F) ^re broad functions of energy. 
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Results 

GCR on a Half- Space 



Figure: Coordinate variables. 


14/25 


355 




Results 

GCR 16 O Fragment flux for the 1977 solar min 


(D 1 = G 1 F 



Figure: The 16 0 fragment flux at (0,0,5). 


16/25 


356 




Results 

GCR on a Circular Cylinder 



Figure: Circular Cylinder 



Results 

GCR 56 Fe Primary flux for the 1977 solar min 


<t>° = G° F 



Figure: The 56 Fe primary ion flux at 

( 0 , 0 , 0 ). 


Figure: The 56 Fe primary ion flux at 

( 0 , 0 , 18 ). 
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Results 

GCR 16 O Fragment flux for the 1977 solar min 


0 1 = G 1 F 



Figure: The 56 Fe primary ion flux at Figure: The 66 Fe primary ion flux at 
(0,0,0). (0,0,18). 
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Results 

GCR 56 Fe Primary flux showing antisymmetry 


<t>° = G° F 



Figure: The 56 Fe primary ion flux at Figure: The 66 Fe primary ion flux at 
(14,0,9) when a!pha=G Deg. (14,0 f 9) when alpha=90 Deg. 
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Results 

GCR 16 O Fragment flux showing antisymmetry 


0 1 = G 1 F 



- D 5 


I « 

? 1 § 
■» a 

, Q .S 


The 16 0 fragment flux at 
(14,0,9) when alpha— 0 Deg. 


Figure: The 113 0 fragment flux at 
(14,0 r 9) when alpha=90 Deg. 
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Conclusions and Future Work 

Conclusions 


0 Voiterra Integral Equation 

0 Solved using a Neumann Series 
solution. 

o Green's functions 
o Closed form approximations 
for 6° and G 1 . 

Showed results for the GCR 
boundary condition. 

0 Future Work 

Figure: Depiction of solar radiation. 

{NASA) 

22/25 



359 





Ac knowled gem en ts 


Acknowledgements 


I would like to acknowledge the NASA support that I received for this work 
under the Graduate Student Research Program (NASA NMXG 9 AJQ 6 H) 
and the Virginia Space Grant Consortium. 




VIRGINIA 
SPACE GRANT 
CONSORTIUM 


23/25 


References 

References 


□ 

B 

□ 

□ 

□ 

B 

□ 

B 


□ 


B 


Committee on the Evaluation of Radiation Shielding for Space Exploration. Managing Space Radiation Risk in the Mew 
Era of Space Exploration, The National Academies- Press, Washington, D.C., 2003. 

Merlens, C, Wilson. J.W. . et a I Coupling of Multiple Coulomb Scattering with Energy Loss and Straggling in 
HZETRN, Advances in Space Research 40. 1357-1357, 2007. 

Miroshnichenko. Leonty I Radiation Hazard in Space. Kluwer Academic Publishers. Boston. MA, 2003. 

Shinn. J.L.. et al. HZE Beam Transport in Multilayered Materials, Radiation Measurements 23, 57-64, 1994. 

Shinn. J.L.. etal. A Green's function method for heavy ion beam transport, Radi at. Environ. Biophys 34, 155-159. 1995. 

Shinn. J.L,„ etal, Multilayer analysis of an Iron radiation beam experiment, NASA TM-4753, National Aeronautics and 
Space Administration, 1997- 

ToWnsend. L-W. Khan, F , and Tripathi, R.K , Optical model analyses of 1 55A GeV argon fragmentation : CrOS5 
section s and momentum distributions. Phys Rev. C 43, 2912-2010, 1903. 

Tripathi. R.K.. Townsend. L.W.. Role of intrinsic width in fragment momentum distributions in heavy ion collisions. 

Phys, Rev C 48, R1775-RI777. 1994. 

Tweed, J , Wilson. J W , Tripathi. R K An improved Green's function for ion beam transport. Advances in Space 
Research 34, 131 1-1318, 2004. 

Tweed, J., et al . Recent Progress in the Development of a Multi-Layer Green's Fucntiort Code for Ion Beam Transport, 
STAIF-SOGS (Space Technology and Applications Forum). Albuquerque, NM, February 10-14. 2003 AIP Conf. Proc. 959, 
993-1000. 2008. 



360 



References 

References 


0 

0 

0 

Q 

0 

□ 

0 

0 

0 


Tweed, J , et al. Validation Studies of the GRNTRN Code for Radiation Transport, ICES 2007-01-3118, SAE 37th 
International Conference on Environmental Systems, Chicago, 2007, 

Walker, S.A. et al.. Validation of the HZETRN code for laboratory exposures with 1 A GeV iron ions in several targets, 
Adv Space Res. 35, 202-207. 2005 

Wilson, J W , Analysis of the theory of high-energy ion transport, NASA TN D-8381, National Aero- nau tics and Space 
Administration, 1077. 

Wilson, J.W., Badavi, F.F. Cucinotta, et a! HZETRN: Description of a Free- Space Ipn and Nucleon Transport and 
Shielding Computer Program NASA TP- 3495, National Aeronautics and Space Administration. 1005, 

Wilson, J W.. Townsend. L.W., Schimmerling, W , et al.. Transport Methods and Interactions for Space Radiations, 
NASA R P-1 257, National Aeronautics and Space Administration, 1991. 

Wilson, J. W... Tweed. J Tai. H. r et al.. A simple model for straggling evaluation, Nucl. Instr. and Meth. in Phys. Res. B 
104. 330-302. 2002- 

Wilson, J W , et al, A closed form solution to HZE propagation, Radiation Research 122, 223-220, 1900. 

Wilson, J.W , et al. A hierarchy of transport approximations for high energy heavy (HZE) ions, NASA TM-41 IB 
National Aeronautics and Space Administration, 1939, 

Wilson, J.W., et al, Mon-perturbative methods in HZE transport, NASA TP-3363, Nationa Aeronautics and Space 
Administration, 1993. 



361 


